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After the discovery of a Standard Model (SM) like Higgs boson, naturalness strongly favors the 
next to the Minimal Supersymmetric SM (NMSSM). In this letter, we point out that the most 
natural NMSSM predicts the following CP-even Higgs H; sector: (A) Hə is the SM-like Higgs boson 
with mass pushed-upward by a lighter Hı with mass overwhelmingly within [ma /2, mn]; (B) 
mu, ~ 2u/sin26 = 300 GeV; (C) H; has a significant coupling to top quarks and can decay to 
Hı Hə with a large branching ratio. Using jet substructure we show that all the three Higgs bosons 
can be discovered via gg > H3 > Hı Hz —> bbévjj at the 14 TeV LHC. Especially, the LEP-LHC 
scenario with Hı ~ 98 GeV has a very good discovery potential. 


PACS numbers: 12.60.Jv, 14.70.Pw, 95.35.+d 


Introduction: Supersymmetry provides the most ele- 
gant solution to the gauge hierarchy problem in the Stan- 
dard Model (SM). In the supersymmetric SMs (SSMs) 
with R-parity, we can not only achieve the gauge cou- 
pling unification, but also have a cold dark matter candi- 
date. Recently, the discovery of a SM like Higgs boson at 
the LHC with mass mp around 126 GeV [1] has deep im- 
plications to the SSMs. Although such a relatively heavy 
Higgs boson mass can be achieved in the Minimal SSM 
(MSSM), it generically incurs a large fine-tuning (For the 
possible solutions, see [2]). By constrast, the next-to-the 
MSSM (NMSSM) with an extra SM singlet Higgs field 
S is strongly favored by naturalness [3], due to origi- 
nally its dynamically solution to the Higgs bilinear mass 
u problem and now the SM-like Higgs boson mass en- 
hancement via the relatively large Higgs trilinear Yukawa 
coupling A in the superpotential and singlet-doublet mix- 
ing effect [4-7]. The natural NMSSM may leave hints at 
the light stop sector, but the search is rather model de- 
pendent [8, 9] and barely has relation with Higgs sector 
(Recent attempt to search for the light stop utilizing the 
properties of the SM-like Higgs boson was done in [10].). 


In the natural NMSSM, the second lightest CP-even 
Higgs boson Hg is indentified as the SM like Higgs boson, 
while the lightest CP-even Higgs boson Hı has dominant 
singlet component. Thus, the Hz mass can be pushed- 
upward via the singlet-doublet mixing effect [4-7]. Such 
a scenario can explain the possible di-photon excess from 
Higgs decays [4, 11, 12] since the significant mixing effect 
reduces the decay width of H> — bb and the light charged 
Higgsino may increase the Higgs decays to diphotons. In- 
terestingly, Hı may be used to interpret the slight LEP 
excess for the Higgs mass around 98 GeV [14] (It receives 
some interest [15, 16] recently.), or the LHC excess for the 
Higgs mass around ~113 GeV [17]. A scenario with two 


light higgs and a low-mass pseudoscalar in NMSSM has 
been discussed in [18]. More noticeable features emerge 
when we take the heavy CP-even Higgs boson H3 into 
account. In this letter, we consider the CP-even Higgs 
sector in the natural NMSSM. We point out that natu- 
ralness implies the H3 mass range my, € [300,600] GeV 
and its significant triple Higgs coupling with Hı and H2. 
Such a Higgs sector structure leads us to investigate the 
discovery potential of the whole CP-even Higgs bosons 
from the process gg > Hs —> H,H2. With jet substruc- 
ture, we show that all three CP-even Higgs bosons H; 
can be probed at the 14 TeV LHC. Our search strategy 
is specially suitable for the LEP-LHC Higgs bosons but 
also applies to the general pushing-upward scenario. 


Light Higgs Bosons in the Pushing-Upward Sce- 
nario: The SM-like Higgs boson can be accommodated 
without recurring severe fine-tuning, and we can show 
that the whole Higgs sector is light. Restricted to the 
Z3—NMSSM, naturalness conditions point to a predic- 
tive parameter space 


A: 0.6— 0.7, tang: 1.3— 3.0, 


= Avs : 100 GeV — 200 GeV, (1) 


where tan £ is the ratio of the vacuum expectation values 
for two Higgs doublets, and « is the singlet cubic coupling 
in the superpotential. Also, « is constrained by perturba- 
tivity, and typically is no more than half of A. The stop 
sector should be sufficiently light, e.g., mz, = mz„ = 500 
GeV, and a flavor safe choice A; = —500 GeV. Their 
concrete values will not qualitatively affect our following 
discussions. 


Importantly, A, can be further determined in the 
pushing-upward mixing scenario. The Higgs mass square 


matrix in the Goldstone basis is 
(M) = M + (m3, — °v?) sin? 28, 
(M3)12 = —5 (mi — vu") sin 48, 
(M3)13 = 5 (M3 sin 28 + 2AKv.) cos 26 
(M3)22 = m% cos? 26 + dv? sin? 28, 


1 
(M5)23 = 5(4A7v — MA sin? 26 — 2x5 sin 28)—, 


s 


2 
(M2)33 = zm sin? 28 (= ) 
v 


b AKU? + KARUs sone sin 26 , (2) 


where M4 = 2dv5(A, + Kvs)/ sin 26 defines the largest 

scale among these elements. Let the orthogonal matrix 
diagonalizing MZ be O: OT Diag(mi,,,m%z,,my,)O = 
M2. The singlet-doublet mixing effect can be approx- 
imately studied by decoupling the entries involving the 
first state. Ref. [4] found that, in the case with a large A 
and small u, the realization of pushing-upward scenario, 
which requires (M?)33 < (M2)22, necessitates a cancela- 
tion to reduce the large non-diagonal element (M2)o3: 


1— (A) /2u+ «/X) sin 28 ~ 0. (3) 


Thus, A) is largely determined by u and tan 8, and to a 
less degree, by x. Then we have 


2u 2 k sin 28 
2 xM? les : 4 
MHs A \ sin 98 x 2 (4) 


Recall that k < A, so, to a good approximation, we get 
ma, ~ Ma ~ 2u/sin26, which is about 2.5, relating 
the H3 mass directly with the weak scale naturalness. 

We now summarize the Higgs spectra in the natural 
NMSSM under consideration. First, all the Higgs fields 
are properly light. Hs and its SU(2), partners, the 
charged Higgs bosons H® and the heavy CP-odd Higgs 
Ag, take roughly degenerate masses Ma. Hə is SM-like 
while H; is even lighter. H; is a SM singlet like and then 
can be allowed by the LEP experiment. Note that my, 
is most likely to fall into the region [mp,/2,m,] with the 
lower bound set by forbidding the decay Hz — H, Ay 
(Ref. [19] considered such case.). Otherwise it tends to 
be the dominant decay mode of Hə. In addition, the 
lightest CP-odd Higgs boson A; also has a mass around 
the weak scale. Moreover, a pair of charginos and three 
neutralinos, consisting of the Higgsinos and singlino, are 
light as well. All of them may be detectable at the LHC 
and here we focus on the CP-even Higgs bosons. 
H,—couplings: The Higgs signals at colliders are sensi- 
tive to their mixing angles whose effects, in a standard 
form, are described by the tree-level Lagrangian: 


VIM% 


M2 biog 
Loves Dri z = Hi ZZ + Tiw AH;,wWT*Ww. 


V2v 


mf = 
I Fy Ff + pije HiH; Ae (5) 


with v = 174 GeV. riv. etc., encode the deviations of 
H; from hsm. For instance, we have 
TRV = Oj2. (6) 


riv = 032, r2,v = O22, 


We also include the triple Higgs couplings, which will 
play a crucial role in the search for Higgs bosons. 

We now present the features of H3 couplings. Firstly, 
note that (M2)j2 is a small entry and we can express it in 
terms of O and m3,,. Since my, is a few times of my, ,, 
then it is not difficult to obtain the upper bound 


O12 = —80, S (M§)12/my, ~ (M§)12/(M3)u, (7) 


where (M2),1 gives the dominant contribution to m p. 
Therefore, the trilinear couplings between H3 and the 
weak gauge bosons are negligibly small. Next, the re- 
duced couplings of H3 to the bottom and top quarks are 
given by 


C3 b = — O;, tan 8 + Oy ~ —O;, tan , 
C34 =O; cot 8 + O12 & O11 cot 2. (8) 


Owing to a relatively small tan £ in the natural NMSSM, 
H3 coupling to the bottom quark is not enhanced while 
its coupling to the top quark is significant. They have 
crucial implications to the collider phenomenology of H3, 
e.g., it can be considerably produced at the LHC by 
virtue of the significant coupling to gluons: 


C3,g =1.03C2,4 — 0.06C2, = O11 cot 2. (9) 


Finally, the triple Higgs coupling H3H2Hj, receives two 
possible large contributions and is given by 


AA, Kk u ) AA) 
PO IY A EE) gg 10 
ane V2 ( A Ay ( ) 


It thus has a large AA) enhancement and leads to H3 > 
Hı Hə decay width at the GeV scale and dominates the 
H, Higgs-to-Higgs decay, as provides the most promising 
discovery prospect for H3 and H4, similarly to Ref. [13]. 

We now turn our attention to the lightest Higgs bo- 
son Hı. Interestingly, the LEP collaboration reported 
(with an signal significance 2.3 ø) a slight excess of events 
for a Higgs boson with mass ~ 95 — 100 GeV [14]. Al 
though our discussions on the Higgs bosons and the en- 
suing search strategy are not restricted to this case, it is 
tempting to interpret Hı as the source of this excess. So 
we have 


Br(H, — bb) 
AE E (01 = 0:25. 11 
cd Brsu(M —> bb) ( ) 


For my < 100 GeV, its decay to bb nearly determines 
its total width. Thus, the LEP requires Cy, vv ~ 0.3 
which is a typical value expected from the mixing Higgs 


sector. 
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FIG. 1: A plot on the o#,-mx, plane, with color code de- 
noting my,. Large inverted triangle points satisfy the LEP- 
LHC scenario. We use NMSSMTools 3.2.1 [20], set 24 cC 


GeV 
(300, 500], Æ$ c [-300, 0], and require [42 C [125, 127] 


GeV 
and signal strengthes R2, ,gg(yy) C [1.4,1.6], Ragg(VV) C 


[1.0, 1.3]. 


Signature and backgrounds: In light of the previ- 
ous analysis, the signature gg > H > H,(— bb)H2(> 
WW?) is promising, where we denote W, as hadronic 
decaying W boson and denote Wọ as leptonic decaying 
W boson. The Wọ will suppress the enormous QCD back- 
grounds. The total cross section is 


(2) Br(H3 — Hı H2) Br(Hı — bb) 
oH, =0.2 T AR a Oe 


0.4 20% 90% 

Br(H2 > WiWn) oar (hsm) 

Ne EAA 12 
28% pb ”” as 


where ¢ = e, u. The numerical results are shown in Fig. 1, 
where a plot of the distribution of o p, on the my, —my, 
plane is presented. It can be seen that its values cluster 
well for a given my, (typically within only a few times), 
in particular for heavier H3. 

We implement the simplified model for Higgs bosons in 
Feynrules [21] to generate the UFO format of the effec- 
tive model for MadGraph5 [22], where the parton-level 
signatures are generated. 

The semi-leptonic tt pair production is the dominant 
background (BG), with the NNLO cross section ~ 240 
pb [23]. The subdominant BG W; + bb+jets has cross 
section depending on the renormalization scale, roughly, 
about 40 pb. Other backgrounds can be neglected in 
our signal region. BGs are generated using MadGraph5. 
To avoid double counting, we adopt the modified version 
of MLM-matching [24] with zqcut = 15 GeV. For the 
latter BG, we include up to 2 additional jets and set the 
k—factor to be 2. 

We use PYTHIA6.420 [25] for decaying particles, 
parton-showering and hadronization. However, in or- 
der to employ the BDRS procedure later, we turn off 
the B—hadron decays in Pythia. The produced objects 
are then converted to the HepMC [26] event format and 
passed to Fastjet 3.0 [27] to cluster the final states. The 
final visible particles are requited to have pr > 0.1 GeV 


and |ņ| < 5.0 which are defined as tracks hereafter. Lep- 
tons from signal events should be isolated, otherwise they 
are combined with the tracks to reconstruct fat jets later. 
Additionally, signal leptons are required to have |7| < 2.5 
and pr > 10 GeV. We take b—tagging efficiency of 70% 
with the other light quark mis-tagged probability 1%. 

We choose the C/A algorithm [29] with radius R=1.4 
and pr>40 GeV to cluster the tracks and form fat jets. 
Following BDRS [28], we first break the hard fat jets into 
subjets 71,2 with masses m,;,,. Next, a significant mass 
drop mj, < umj with u=0.667 and not too asymmet- 
ric splitting, i.e., y =min(p7 j» Pij ARS, p/m > Yout 
with Your=0.09 (AR, jẹ is the angular distance), are re- 
quired. If the above criterion are not satisfied, we will set 
j = jı and go back to decomposition. Finally, we filter 
the Higgs neighbourhood, resolving the fat jets on a finer 
angular scale Rai, = min(0.35, R;,;,/2) and taking the 
three hardest objects, with the remains identified as the 
underlying events contamination and hence dropped. 
Events selection and results: Two basic cuts are im- 
posed to trigger our events. Firstly, at least two filtered 
fat jets are required. One of them has two leading subjets 
which pass b—tagging and satisfy |n| < 2.5, and then is 
identified as the H,—jet. Among the remaining fat jets, 
the one with highest pr is regarded as the W;,—jet [30]. 
Secondly, the events must contain exactly one isolated 
lepton. 


NM, 


FIG. 2: Distribution for trigged signal and background of 
PTH, MH,, Mjjtv, AQij, ARg op, Mc. The number of events 
have been normalised to 14 TeV 500 fb~*, and the signal is 
400 times amplified. 


For illumination, here we will take a benchmark point 
inspired by the LEP-LHC Higgs scenario: my, = 98GeV, 


MH, = 125 GeV as well as my, = 400 GeV. Figure 2 
shows the distributions of some important kinematic vari- 
ables. In terms of the plots, we display the cut flow: 


e Cutl: The relatively large mass splitting between 
H; and Hı gives Hı a boost. Therefore, we require 
Prop > 150 GeV, prjjev > 120 GeV, and [PT = 
PT jjev| < 20 GeV. 


e Cut2: It is observed that the longitudinal momen- 
tum of the neutrino from W decay is generically 
small, and hence my,,, can be approximately re- 
constructed by assuming pz, = 0. Practically, cuts 
based on this assumption are sufficiently good. So 
we impose: 95 GeV < my, < 100 GeV, mjje, < 
150 GeV, and mp5; 50, < 440 GeV. 


e Cut3: Because Hə has spin-0 and W only couples 
to the left-handed fermions, the lepton from Wọ 
will align with one of the jets from Wp, decay. It 
allows us to impose a cut |Agy;| < 1.5, namely 
the azimuthal angles difference between the signal 
lepton and (one) jet being sufficiently small. 


e Cut4: The filtered Hy-jet actually contains three 
subjets, the bb and a radiated gluon. So the H-jet 
and its bb subsystem must have a very small angle 
distance. By contrast, the angle distance between 
the Hı-jet and W)-jet is much larger. Thus, we 
require ARņp, 45 < 0.01, and 2.6 < ARy,w, < 3.4. 


e Cut5: We also impose the cluster transverse 
mass of decay product of the He: Mc = 


Phase + Mja + Br < 220 GeV. 


With the above cuts, we obtain the signal significance 
4.42 o excess for the LEP-LHC benchmark point at 14 
TeV 500 fb-!. The cut efficiency and the signals are 
presented in Table I. 


tt W (— lvjj)bb + jets 


Signal 


Total | 1.2 x 10° 1.91 x 107 1.25 x 104 
Trigged]4.95 x 10° 1.45 x 10° 1456.75 
Cutl [3.77 x 10° 1.61 x 10° 639.5 
Cut2 1932 203 119.75 
Cut3 1512 155.2 105.5 
Cut4 108 47.75 56.25 

Cut5 84 47.75 55 


TABLE I: Number of events after each cut for background 
and signal (normalized to 500 fb~'). The signal significance 
S//S + B has reached to 4.02 and with the precise 4.42 ø 
excess for the LEP-LHC benchmark point. 


Since Fig. 1 shows obvious cluster behavior, the whole 
parameter space with pushing-upward effect can be ex- 
plored. Using BDT analysis [31], we consider six repre- 
sentative points to demonstrate the search prospect, and 


the discovery signal significance for each case is given 
in Table II. Some observations can be made: (A) For a 
given mz, a lighter Hı shows better discovery potential; 
(B) Increasing H3 mass helps to boost H, but the cross 
section is reduced. Thus, a moderately heavy H3 ~ 400 
GeV and relatively light Hı have the most promising dis- 
covery potential; (C) Most of the parameter space is dis- 
coverable except for simultaneously light H3 and heavy 
A, e.g., the benchmark point B1, despite of a rather 
large cross section, has a quite low signal significance. 
The situation can be further improved when we take 
H= into account. H~* can be produced associated with 
a single top, with a moderately large cross section at the 
small tan 8 region. Moreover, it can decay to Hı and 
W with a substantial branching ratio and hence provide 
a way to probe Hı and H=. In this case H* can be 
produced with a larger pr and the boost will be easier for 
lighter H~, not very sensitive to my,. So it can provide 
a complementary or even more promising channel for the 
pushing-upward scenario. We leave it for the future work. 


mu, (GeV)|mu,(GeV)|}o (fb) | -= 


VS+B 
B1 100 300 70 | 0.81 
B2 65 300 50 | 3.84 
B3 98 400 25 | 4.73 
B4 65 400 20 | 7.68 
B5 100 600 2 2.79 
B6 65 600 2 4.99 


TABLE II: Discovery signal significances for 6 representative 
points at 14 TeV 500 fb™t. We design 25 kinematic variables 
for BDT analysis [31]: Ér, pr,w, Mw, Njet, PT,bı, PT.by PTL; 
MT ev, PT,wjt, PT,wjr, PT jjev; ARej, Adiw, PT.H3, Mev, Ew, 
and mH3- 


Conclusion: We pointed out the specific features in the 
CP-even Higgs sector of the natural NMSSM, and showed 
that all three CP-even Higgs boson H; can be probed at 
the 14 TeV LHC. 
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